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Based on the properties of homogeneous functions of first degree, we derive in a mathematically 
consistent manner the explicit expressions of the chemical potential /i and the Clausius entropy S 
for the case of monoatomic ideal gases in open systems within phenomenological thermodynam- 
ics. Neither information theoretic nor quantum mechanical statistical concepts are invoked in this 
derivation. Considering a specific expression/ value of the constant term of S, the derived entropy 
coincides with the Sackur- Tetrode entropy in the thermodynamic limit. We demonstrate however, 
that the former limit is not contained in the classical thermodynamic relations, implying that the 
two aforementioned entropies are distinct in finite systems. 
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I. INTRODUCTION 



The first law of thermodynamics presents an expression of the energy conservation principle, 
differentials the former law is given as 

TdS = dU + PdV - fidN , 



In terms of total 



(1) 



where T, S, U, P and V denote the temperature, the (Clausius definition of) entropy, the internal energy, the pressure 
and the volume, respectively. N denotes the number of particles and /x the chemical potential associated with the 
particle exchange (for one species of particles). A possible solution of Eq. ([T]) for the case of monoatomic ideal gases, 
with respect to the entropy as a function of U, V and N, has been initially computed within quantum statistical 
mechanics, with a Hamiltonian consisting only of kinetic energy contributions, presenting a solution of Boltzmann- 
Gibbs statistics as well. The obtained result is known as the Sackur- Tetrode (ST) entropy and it has the following 
expression 



S ST (U,V,N)=limJ B ln( J ^ 



2 \4TTm J 



(2) 



where k B and h are the Boltzmann and Planck constants respectively, and $ is the 67V-dimcnsional phase-space 
volume of the region enclosed by the constant energy hyper surface The term m denotes the mass of the 
particles. The factor N\ stems from the traditional way of resolving the Gibbs paradox interpreted by means of 
the quantum mechanical concept of the particles indistinguishability, taking in this way into account the invariance 
under permutations. An information theoretical derivation of the ST-entropy has been also demonstrated by A. 
Ben-Nairn in his textbook where the former function is the result of the combination of the following four terms: 
positional uncertainty, momenta uncertainty, quantum mechanical uncertainty principle and the indistinguishability 
of the particles. In the same book one may also find an alternative derivation of the factor TV!, with the Schrodinger 
equation being the point of departure. Sst is then expressed in terms of the de Broglie wavelength. 

From the mathematical point of view, the ST-entropy in Eq. ([2J presents a homogeneous function of first degree. 
We remind that a function / : M. q + — > R of class C n - 2 is said to be homogeneous of rth degree (rd) generally, when 
it satisfies the following relation 



/(ax) = a T /(x) 



(3) 



where a G M+ and x G R+. Then, homogeneity of degree one (Id) is characterized by Eq. © for r = 1. Requiring 
the aforementioned mathematical property as a fundamental property of the thermodynamic entropy S, we aim in 
this paper to derive a solution of Eq. ([l} for monoatomic ideal gases in open systems within the context of the 
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phcnomcnological thermodynamics. It is demonstrated that the entropy function S obtained from our homogeneity 
analysis coincides with the ST-cntropy for a properly chosen constant term of S. In the current derivation however, 
we do not explicitly invoke the N — > oo limit, in contrast to the Sst derivation, as can be seen in Eq. (0. We further 
show that the first law of thermodynamics does not contain the thermodynamic limit, implying that if ld-homogeneity 
is a fundamental entropic property, then the monoatomic ideal gas entropy is not given by S ST in Eq. ([2]) but is given 
by the Clausius entropy S. This result is of great theoretical interest, since it indicates that, under the request of the 
entropic ld-homogcncity, the statistical background of S ST does not succeed to derive the correct entropy expression 
of monoatomic ideal gases. 

The paper is organized as follows. Section [TT] reviews the basic relations characterizing the homogenous functions of 
first degree. After having outlined the formalism of the former functions, in Section Unl we apply it to thermodynamic 
functions deriving the chemical potential and the entropy expression for monoatomic ideal gases in open systems. 
Conclusions are presented in Section IIVI 



II. HOMOGENEOUS FUNCTIONS OF FIRST DEGREE 



Let us consider a Id- homogeneous (ldH) function / : M. q + — > R of class C"- 2 . Then, according to Eq. (|3]) (for 
t = 1), /(x) satisfies the relation 

/(ax) = a/(x). (4) 

An example of the former function for the case q = 2 can be provided as f(x,y) = xexp(x/y), where x% = x and 
X2 = y. The differentiation of Eq. (HJ with respect to a yields 

Z^i-H? T = /( x - (5) 

Eq. ([5]) is valid for any a £ R+, thus we can choose a = 1. We then obtain 

/(x) =5> 4 A,(x), (6) 

i=l 

where Aj(x) := df(x)/dxi. Regarding our example above, we have -Ai(x) = A x (x,y) = (1 + x/y) exp(x/y) and 
y4 2 (x) = A y (x,y) = — cxp(x/y)x 2 /y 2 . Indeed then, f{x,y) can be expressed as in Eq. ©, i.e., f(x,y) = xA x {x,y) + 
y A y (x,y). The aforementioned equation is referred to as Euler's ld-homogeneous function theorem, stating that a 
ldH- function /(x) can always be expressed in the form of Eq. It can easily be shown that the relation between 
Eq. ([4]) and Eq. ([6]) is bijective. From the latter equation we see that the functions ^Lj(x) are Od- homogeneous (OdH), 
i.e., Ai(ax) = >lj(x), which is verified also from our example, A x (ax,ay) = A x (x,y) and A y (ax, ay) = A y (x,y). 
Computing the partial derivative of /(x) in Eq. with respect to Xj, we obtain 

»tV-W + Y.^, (T) 



dxj 4—1 dxj 

j / = i 



implying 



■ 

Xi - ■ • =0. (8) 



x ^ 8Ai(x) 



4-1 " 9x 



As can be seen, Eq. ([5} is satisfied by any ldH-function and thus it presents another characteristic relation of the 
former class of functions. In fact, Eq. ([8} presents a set of q conditions, since j = 1, . . . , q. Due to /(x, y), the former 
conditions take the form 

dA x (x,y) dAJx.y) fx 2\ f x\ ( x 2 2x\ fx 

* , v ; + y -k = x — + - cx p - + v - — - — ex p - 

ox Ox \y z yj \yj \ y A if) \y 

dA x (x,y) dAJx,y) ( x 2 2x\ ( x\ fx 3 2x 2 \ fx\ 

' y V -5 =4 -Tff-Ta exp - )+y(— + — exp - =0. (10) 



oy oy \ y A y J \vj \y y A J \y 



3 



Now, we want to prove the reverse, namely if Eq. ((8j) is true, then /(x) is ldH. Eq. ((8j) is valid for any translation 
/(x) of /(x), i.e., /(x) = /(x) + c, where c is a constant, so that Aj(x) := df(x)/dxi = df(x.)/dxi = A 4 (x). 
Multiplying Eq. © by dxj and then summing over all j's in the aforementioned equation, yields 



oxj 1-— 1 ox 

J = l 2—1 J 2 — 1 J — l 



The change in the order of the summation in Eq. (jTTJ) can be performed, since 2j=i(Xi=i x i [dAi(x)/dxj\ dxj = < 
oo) < oo. Adding the term Xh=i Aj(x)c&Ki in Eq. (TTT]), we obtain 



^Aj(x)efccj =^A i (x)da; j +^a; j cL4 i (x) d/(x) = rff ^^A^x) 

2=1 2—1 2—1 ^ 2 — 1 ' fl2) 

/(x) = ^a^-A^x) + a => f(x) = ^2xiAi{x) + (a-c), 



where a is an integration constant. We can always choose c = a, so that Eq. (fT2)) is identified with Eq. ©. Herewith, 
we have proven that the conditions in Eq. ([8]) are necessary and sufficient conditions for having ld-homogeneity. 
Thus, they can be considered as an alternative definition of ldH-functions. An important consequence of Eq. (JSJ) (or 
equivalently Eq. (TTT]) ) is that the differential of the function /(x) is exact 

2—1 2—1 2 — 1 2 — 1 1 

implying simply, that there are indeed functions satisfying Eq. |QJ. 

It is worth remarking that Eq. (J8j) , because of the equality dAi(x.)/dxj ~ dAj(x)/dxi which holds true for any 
C"- 2 -function, can be rewritten as x i \dAj{x) / dxt\ = 0. One can easily verify that the former equation is the 
respective result of Eq. $Q for OdH functions. It becomes obvious then that Eq. © fully characterizes both the 
ldH-function /(x) and the OdH-functions Aj(x) simultaneously. 

III. MONOATOMIC IDEAL GAS ENTROPY OF THREE INDEPENDENT VARIABLES 

We now make use of the results presented in the previous section on ldH- and OdH-functions to model thermo- 
dynamic quantities. Considering three independent variables, i.e, U , V and N, representing the internal energy, the 
volume and the number of the particles of a system under scrutiny, respectively. Further, we assume a ldH-function 
S = S(U, V, N) as the entropy. Then, according to Eq. ([5]), S can be expressed as follows 

b[U > V > N) - T(U,V,N) U+ T(U,V,N) V T{U,V,N) N ' (14) 

where 

1 dS(U,V,N) P(U,V,N) dS(U,V,N) n(U,V,N) dS(U,V,N) 



T(U,V 7 N) ' dU T{U,V 7 N) ' dV T(U,V,N) ' ON 



(15) 

In thermodynamics, Eq. (|14[) is referred to as the Euler equation. The left hand side of the definitions in Eq. (|15[) arc 
concrete physical quantities. T, P and \i denote the temperature, pressure and chemical potential, respectively. They 
suppose to be OdH-functions, otherwise the assumption of S being ldH cannot hold true. As shown in the previous 
section, S has to satisfy the conditions in Eq. (|5]). In the thermodynamic notation the former conditions take the 
form 

"w^w^U-"- < 16 *> 

d 1 d P da 
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If the OdH- functions P,T and a satisfy Eq. ((16)) , then the differential of S is exact yielding the hrst law of thermo- 
dynamics, given in Eq. ([T]). Considering the relations in Eq. (|16|) simultaneously, they can be written as 



Ud 



Vd 



Nd 



0. 



(17) 



We identify Eq. (|17[) with the well known Gibbs-Duhem relation 

Let us now explore the conditions in Eq. (|16[) for the case of an ideal monoatomic gas in an open system, dN ^ 0. 
The expressions of pressure and temperature are given by 



P(U,V,N)=P(U,V) 



2U 
W 



T(U, V, N) = T(U, N) 



2U 
3fc R A 



(18) 



The pressure in Eq. (fT5|) has been derived within the kinetic gas theory. The temperature in Eq. (jTHJ) , within classical 
thermodynamics, presents an ad hoc definition, being compatible with the equation of state for the ideal gases, i.e., 
PV = k B NT. Indeed, combining P and T in Eq. (fT8|) one obtains the former equation of state. Apparently, the 
functions P and T are OdH. In the mathematically trivial case of constant pressure and temperature, dP = = dT => 
U ~ N and V ~ N, Eq. (|16[) yields a constant chemical potential /j,. In this case, the respective entropy function in 
Eq. (Q3J) is ldH and linear with respect to N, i.e., S(U,V,N) = S(N) <~ N. For the general case of varying pressure 
and temperature i.e. dP ^ and dT ^ 0, from Eq. (fT6)) . we read 



±,,(U,V,N) = ^(U,V,N)-±, 
d , . 2U 

w« u > v >w = -wn> 

f) 1 F>rr 

wm v,N) = -- m v,N ) + — 2 , 



(19a) 
(19b) 
(19c) 



respectively. Solving this system of differential equations, we arc able to determine the expression of the chemical 
potential compatible with Eqs. (JXJ) and (fT5|) . The aforementioned system can be solved stepwise. We first solve Eq. 
(|19al) . The result will depend on an integration constant Z(V, N), with respect to U . Then, plugging the former result 
into Eq. (|19b[) . Z(V,N) reduces to an integration constant Z'(N). Finally, following the same procedure with Eq. 
(I19c[) , Z'(N) reduces to a constant, which cannot be further reduced. Thus, for varying pressure and temperature, 
the monoatomic ideal gas chemical potential is determined apart from a constant £, having the expression 



H(U,V,N) 



u_ 

N 



2 (N 



As expected, the ^-function is OdH. The respective entropy function S in Eq. (|T4")) is then ldH, given as 



S{U,V,N) = k B N 









5 3 ; 


Hi) 






+ 2-2^. 









(20) 



(21) 



Studying the validity of S, we observe that the function S in Eq. (|21j) can be derived by means of the ld-homogeneity 
formalism when P,T, /i ^ {0, const., ±oo}, beforehand, and when dU ^ 0, dV ^ and dN ^ 0, simultaneously. If one 
of the former differentials is equal to zero, implying that the respective variable is a constant, then the ld-homogeneity 
formalism is not applicable, since in this case not all of the partial derivatives of S, i.e., 1/T, P/T or fi/T arc OdH. 
Comparing now the Clausius ideal gas entropy S in Eq. (f2Tj) with the Sackur- Tetrode entropy S ST in Eq. ([2]), we 



observe that they both have the same form when £ takes the value £ = In (J^— )• This specific expression of £ cannot 
be obtained by virtue of the ld-homogeneity analysis but only within quantum statistical mechanics. However, even 
then, it is merely a constant and has no physical impact on the entropy expression. So, we observe that the physical 
part of Sackur- Tetrode entropy coincides with the one of the ideal gas entropy S in the N — > oo limit. We stress 
that the analytical expression of ST-entropy is obtained after the explicit computation of the former limit, as can be 
seen in Eq. ([2]), while S has been derived in practice for any N > 0. Thus, the two aforementioned entropies may 
be identified only under the assumption that the thermodynamic limit is implicitly included in the thermodynamic 
relations. This in turn would imply taking into account the restrictions mentioned above regarding the OdH-functions 
T, P and /i, that as much the energy U as the volume V are considered in the same limit. If this is the case, then one 
has to be very careful in discussing phcnomcnological thermodynamics of finite systems. In different case, studying 
ideal gases in open systems for finite ./V we have S ^ S ST , yielding merely the same picture in the — > oo limit. 
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A possible way to shed light on the aforementioned issue, namely whether phenomenological thermodynamics 
implicitly contains the thermodynamic limit, or equivalently whether S and S ST are identical, is by considering the 
general solution (gs) of the differential FLT in Eq. ([I]), regarding the ideal gas entropy function. The latter is given 
as in Eq. (j2"Tj) , with the £-term being a function of N, i.e., £ — > Q(N), namely 



Sgs{U, V, N) = k s N 



1 N / 2 V A 



---6(A) 
2 2 v ; 



(22) 



In this general case, the chemical potential, fj, gs = —T(dS gs /dN), may generally not preserve the OdH property. From 
Eqs. © and (f2"2"j). we observe that in the case of ST-entropy the 0-function in finite systems is inhomogeneous and 
has the explicit expression 



e sT (A) = £ + ■= 1 - ln(JV) + 



ln(AT!) 



A 



(23) 



satisfying the limjv->oo Qst(N) = £ limit. Let us then explore which is the most general expression of 0(A) being 
compatible with ld-homogeneity, S gs — > (and /j, gs — > ^gg H ). Therefore, we compute again the conditions in Eq. 

(|16p. in terms of S gs now. From the first two conditions we do not obtain any new information, yet the last one yields 



A 



<9 2 0(A) d®(N) 







0(A) = a t 



02 

A 



dN 2 dN 

where a\ and a,2 are integration constants. Plugging the result in Eq. (f24|) into Eq. (f22|) we obtain 

1 „ P(U,V,N)„ ^(U^N) 



(24) 



c-ldH 



{U VN) = - u+ WW y 

[ ' ' ' T(U,V,N) T(U,V,N) 



T(U,V,N) 



-TV- 



3, 



k B N 



hi 



In 



(25) 



-k B a.2 



where /i°^ H = (i in Eq. ([20]) with ai = £. 



Comparing Eqs. (|2ip and (|25[) . it becomes obvious that S g ^ K contains an 
extra additive constant. However, in Eq. (|12j) we showed that such a constant can always be transformed away. Thus, 
the |fcBa2-term has neither physical nor mathematical impact on so that the former is identified with S in Eq. 
(|2ip . Eq. (|24[) then, unveils that the thermodynamic limit is not included in the phenomenological thermodynamic 
equations. If it would be included, then the a2-term would not appear, limjv-xx, (12/N = 0. Accordingly, the Clausius 
IdH entropy function S in Eq. (|2ip . derived by means of purely thermodynamic relations, is distinct from the IdH 
ST-entropy, derived within statistical mechanics in Eq. ([2]). 



IV. CONCLUSIONS 



Studying homogenous functions of first degree (Id), we demonstrated that they can be defined through a set of 
g-specific conditions given in Eq. ([5J , where q denotes the number of independent variables of the aforementioned func- 
tions. If the ld-homogeneous function under consideration is the thermodynamic entropy, then the former conditions, 
expressed compactly in a single equation, are identified with the well-known Gibbs-Duhem relation. 

Considering particularly the pressure P and the temperature T for the case of monoatomic ideal gases in open 
systems, where the number of particles N, the internal energy U and the volume V are the independent variables, 
we derived by means of the ld-homogcncity conditions, the analytical formulas of the associated chemical potential 
/i, related to the particle exchange, and the associated thermodynamical entropy S. Interestingly enough, for dP ^ 
and dT 7^ (=> dfi ^ 0), the former entropy coincides with the well known Sackur- Tetrode entropy S ST , apart 
from a non-physical A-translation. In practice however, distinct physical conditions are required for their derivation. 
In fact, the two aforementioned entropies can be identified to each other only under the assumption of a "hidden" 
thermodynamic limit in the phenomenological thermodynamic relations. We explicitly showed yet, that the former 
assumption does not hold true. Accordingly, the Clausius entropy S, derived on the ground of homogeneous functions 
of degree one within phenomenological thermodynamics differs from the Sackur- Tetrode entropy S ST , derived within 
statistical mechanics, sharing merely the same expressions in the thermodynamic limit. 
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